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In order to understand how screening is modified by elec- 
tronic interferences in a mesoscopic isolated system, we have 
computed both analytically and numerically the average ther- 
modynamic and time dependent polarisabilities of two dimen- 
sional mesoscopic samples in the presence of an Aharonov- 
Bohm flux. Two geometries have been considered: rings and 
squares. Mesoscopic correction to screening are taken into 
account in a self consistent way, using the response function 
formalism. 

The role of the statistical ensemble (canonical and grand 
canonical), disorder and frequency have been investigated. 
We have also computed flrst order corrections to the polaris- 
ability due to electron-electron interactions. 

Our main results concern the diffusive regime. In the 
canonical ensemble, there is no flux dependence polarisability 
when the frequency is smaller than the level spacing. On the 
other hand, in the grand canonical ensemble for frequencies 
larger than the mean broadening of the energy levels (but still 
small compared to the level spacing), the polarisability oscil- 
lates with flux, with the periodicity h/2e. The order of mag- 
nitude of the effect is given by 5a/ a oc (Xs/Wg), where A is 
the Thomas Fermi screening length, W the width of the rings 
or the size of the squares and g their average dimensionless 
conductance. This magnetopolarisability of Aharonov-Bohm 
rings has been recently measured experimentally M and is in 
good agreement with our grand canonical result. A prelimi- 
nary account of this work was given in H . 



I. INTRODUCTION 
A. Motivation of the paper 

Transport and thermodynamic properties of meso- 
scopic metallic samples are known for a long time to be 
quite sensitive to the quantum phase coherence of the 
electronic wave functions at low temperature 0] . In par- 
ticular, it was shown that a mesoscopic ring threaded by a 
magnetic field exhibits a persistent current periodic with 
the flux and with a period of one flux quantum (/)o = h/e. 
This current was measured in several experiments [^ |^] 
H, corresponding to different situations and materials. 
More recently, the response of a mesoscopic system to a 
time dependent magnetic field was studied both theoret- 
ically and experimentally [^. 

The purpose of the present paper is to investigate how 
screening is infiuenced by electronic interferences in a 
phase coherent sample. In this spirit, we discuss how the 
electrical polarisability a, i.e. the response of a metallic 
sample to an electrostatic field, a quantity which is di- 
rectly related to screening, is also sensitive to mesoscopic 



effects, a can be experimentally measured by coupling 
the mesoscopic samples to a capacitor. For small volt- 
ages, corresponding to the linear regime, the change in 
the capacitance is directly proportional to the average 
polarisability < a > of the particles. 

When the typical size L of the particle is much larger 
than the Thomas Fermi screening length As, a is mostly 
determined by the geometrical shape of the particle. For 
a spherical particle, it is proportional to the volume V of 
the particle, with a small correction of the order of Xs/L. 

Now, what happens in the case of a quantum coherent 
sample, (i.e when L is the order of the phase coherence 
length)? Is there a mesoscopic correction to the polaris- 
ability, due to contribution from electronic interferences? 

In order to answer this question, we have calculated 
the flux dependent correction to the average polaris- 
ability, which we call "mesoscopic correction" , for two 
dimensional mesoscopic samples in the presence of an 
Aharonov-Bohm flux. Two different geometries have 
been considered, namely a two dimensional ring or square 
in an in-plane electric field. The role of statistical en- 
semble, disorder, and finite frequency have been investi- 
gated. In the case of a time dependent electric field, we 
distinguish two different limits, whether the frequency 
is smaller or larger than the inverse relaxation time, or 
mean level broadening 7. The case w ^ 7 corresponds 
to the thermodynamic limit whereas w ^ 7 corresponds 
to the finite frequency limit. 



B. Scope of the paper 

This paper is organized as follows: the definition of 
the electrical polarisability and the classical result for a 
macroscopic sample are given in section O. 

In section [II, we derive the expression of the electri- 
cal polarisability of a system of non-interacting electrons, 
in the spirit of the early work of Gorkov and Eliashberg 
(GE) M. A giant unphysical polarisability is obtained. 
However it is possible to discuss in this very simple model 
the flux dependence of two quantities which contribute 
to the polarisability, namely the matrix element of the 
position (or screened potential) operator and the contri- 
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bution of energy levels. Then, we examine in section 
how screening can be taken into account. In particular 
we show how to incorporate contributions of electronic 
interferences. Using those results, a general expression 
for the polarisability and its mesoscopic correction for a 
time varying electric field are derived in section |V|. 



We study then separately the case of the grand canoni- 
cal and canonical statistical ensemble. The grand canon- 
ical average polarisability is found to depend only on the 
flux dependent matrix elements of the screened poten- 
tial whereas the Canonical average depends also on the 
energy levels. As a result when the frequency is smaller 
than the level spacing the magneto-polarisability is found 
much smaller in the canonical ensemble compared to the 
grand canonical case. Finally, we present results on the 
Hartree-Fock correction to the zero frequency polarisabil- 
ity. Section VI is devoted to a comparison of these results 
with recent experiments. 



II. BASIC CONCEPTS 
A. Screening and polarisability 

When an external field Eext = — gi"ad((/iext) is apphed 
on a conductor, the charge density adjusts in such a way 
that the effective field inside the conductor cancels. This 
phenomenon is called screening. In principle the compu- 
tation of the potential inside the sample requires to solve 
a complex N body hamiltonian taking into account the 
long range Coulomb interaction between electrons. How- 
ever it is possible in general to treat the interactions in a 
mean field approximation. The screened potential inside 
the conductor is then given by: 



ri 



t(ri)+ / (5/9(r2)f/(ri,r2)dr2 



(1) 



Where 5p{r2) is the charge density induced at the point 
r2 and C/(ri,r2) is the coulomb interaction potential. 
In the linear response limit, Sp{r2) can be described by 
the electric response function: x(ri,r2)7 which relates 
the variation of charge in the system at r2 to a local 
perturbation at ri:(5(ri). 



<5/o(r2) = / x(ri,r2)0(ri)dri 



(2) 



Relation m can then be written more conveniently in a 
matrix form: 



= (i-c/x)- 



(3) 



In a clean infinite system, the induced charge density 
is simply related to the screened potential (j) in Fourier 
space: Sp{q) = x{<l)4'{<i)i where <j){q) = J e"''''(/)(r)dr. 
The problem is more complex in a disordered system, 
which does not have the translational invariance. In any 
case for a finite system, it is possible to compute the 
response function from the eigenstates ja > of the system 



in the absence of the external field, 
eigenvalues ta- 



associated to the 
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a=l I3<a 



C(i"i)V'/3(r2)V'a(ri)V'^(r2) 



e/3 



-I- c.c. 



(4) 



In most cases we will identify these eigenstates as the so- 
lution of the hamiltonian of the system without electron 
interactions which will be taken into account only by con- 
sidering the screening on the applied field. The effect of 
interactions on the eigenstates will be however discussed 
in section VII within the Hartree Fock approximation. 

The application of an electric field on a finite metallic 
system, results in an induced dipolar moment P which, 
for a sufficiently small field R, is proportional to the ap- 
plied field. 



P = / p (r) rdr = aE 



(5) 



where a is by definition the polarisability of the conduc- 
tor, which can be expressed as a function of the response 
function: 



a = -^ / x(ri,r2)E.r2^(ri)dridr2 (6) 



B. Classical polarisability 

1. Plate with perpendicular electric field 

The polarisability of a macroscopic piece of metal is 
described by classical electrostatics. In this approxima- 
tion, the dipolar moment can be calculated by solving 
Poisson equation, AV -f — = 0, in all space and impos- 
ing as well that the conductor is at a constant potential. 
The problem is solvable analytically for a sphere of vol- 
ume V: a — eo{4:Tr/3)V. It is also possible to estimate 
the polarisability of a metallic plate of thickness a in a 
constant electric field perpendicular to the plate. The in- 
duced charge density is such that the electric field cancels 
inside the plate. Neglecting border effects, we have: 



Eeo 



'(')^^K^-i)-^(-i) 



(7) 



Integrating equation pi one finds the classical Polaris- 
ability: 



a = eoa 



eoV 



(8) 



^This point will be discussed later on 
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eoL^ln{D/L) [^. Similarly as in the 2d case, the charge 
density has to spread over the entire wire in order to 
screen the external potential. 

The case of a ID ring is also exactly solvable. In order 
to take screening into account, we first assume that the 
induced charge density in the ring is of the form p{9) = 
Xcos{9)6{z)S{p — R), from which we deduce the induced 
potential in the ring: 
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p{e') 



Stteo sin I 



rde' 



(10) 



The integral is logarythmically diverging and it is nec- 
essary to introduce a cutoff Be — W/R, related to the 
finite width W of the ring: (Pm{0) « _A£2!MMM. 

The screened potential reads (j){9) — cos{9) [E + AJ), 

with J = j^ In ( ^ J . The polarisability is obtained 

either from the potential calculated above, either from 
the expression of the charge. This self-consistent relation 
sets the value of the parameter A, from which we deduce 
the polarisability of the ring. 



FIG. 1. Induced charge density of a conducting strip in an 
in-plane uniform electric field: the function diverges logary- 
thmically as a function of the cutoff A/a. 
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2Jao 



(11) 



where ao = cqit R which reads in the limit J ^ 1: 



2. 2d and Id systems in longitudinal electric fields 

We now examine the case of 2d conductors, and more 
specifically a strip of length L and width W, with L 3> 
VF, in an in-plane electric field perpendicular to its long 
axis. In order to determine the induced charge density, 
the effective electric field inside the strip is imposed to 
be zero. In the limit of an infinite strip, this condition 
leads to the following equation: 



eoTT 
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(12) 



As a result both longitudinal polarisabilities of a ring 
and a disk scale like the cube of their radius, just like the 
polarisability of a sphere. 



C. A mesoscopic correction to the polarisability ? 



1 PJ^') 
27reo X — x' 



dx' 







(9) 



We have calculated numerically the solution of in 
the case of an infinite strip (fig. nh. In contrast with 
the 3 dimensional case, the charge distribution spreads 
over the full width of the strip in order to cancel the 
external field. The numerical solution is very close to 
a logarithm. Assuming that the density is of the form 
p{x) « EW log (-f - x) / log (-f -I- x) , the integral | of 
the previous function gives a ~ egW^L. 

This result extrapolates to a 2 dimensional sample of 
typical size W, for which the in-plane polarisability scales 
like a ~ eoW^. For instance the polarisability of a disk 
of radius R is a = eofR^ 0. 

In the case of a one dimensional wire, whose length L 
is much longer than its diameter D, in an electric field 
along its axis, the polarisability has been shown to be 



In the following, we define the quantity we wish to 
calculate. As we have just seen, the classical polaris- 
ability of a metallic sample is determined by geometri- 
cal factors. However, due to the Pauli principle, charges 
can not strictly accumulate on the border of the sample, 
but the charge density distribution rather extends on the 
screening length A^. This effect reduces the polarisabil- 
ity compared to its classical value by a quantity SaxF 
of the order of Xs/L, where L is the typical size of the 
conductor along the electric field. In a coherent sample, 
electronic interferences might also give rise to a contri- 
bution Sag to the polarisability. The total polarisability 
can thus be written: a — ad + SaxF + ^ckq. In order to 
show evidence of this correction, a magnetic fiux can be 
applied through the sample in order to introduce a flux 
dependent part in the phase difference of electronic tra- 
jectories. We will call "mesoscopic correction" the fiux 
dependent part of the polarisability Sa{^). We will focus 



on the ensemble average of this quantity (5a($), which is 
$0/2 periodic. We will characterize the flux dependence 
of la by d^{a) = ^($o/4) - ^(0). 



D. Model used for the numerical simulations 

Part of the analysis we have done rely on computer 
simulations. The two dimensional rings or squares are 
modelled by the Anderson model. Each atomic site is 
coupled to its nearest neighbour by a hopping term t. 
Disorder is introduced by on-site energies randomly dis- 
tributed in the interval [—w,w]. The Anderson hamilto- 
nian can be written in second quantization: 



H 



J2 ^^'^l'^^ 

k 



te"f"'clck+i 



''^ '■^k+l'^k 



(13) 



where cj. is the creation operator associated with site 
k. The magnetic flux is introduced in this hamiltonian 
through the phase factor of the hoping element with 
(j)k — 27r$/$o {xk — Xk+i)/L for a ring of size L. In 
the case of squares we have taken into account the pen- 
etration of the magnetic field in the sample: the phase 
factor appearing in the hopping matrix elements is com- 
puted from the integral of the potential vector between 
concerned neighbors. In two dimensions, an estimation 
of the elastic mean free path is given by le « 30 (;^) Q] 
allowing us to choose values of t/w so that the system is 
in the diffusive regime. In order to average quantities we 
take advantage of ergodicity properties in the diffusive 
regime and average over the number of electrons in the 
system between one quarter and three quarter filling. 



III. QUANTUM POLARISABILITY OF A 
SYSTEM OF NON-INTERACTING ELECTRONS 

We now examine the electrical response of a quantum 
system of non- interacting electrons, which was first cal- 
culated by GE [gj for a metallic grain both in the diffusive 
and ballistic limit. Their conclusion is quite surprising: 
in the diffusive regime, a is bigger than the classical po- 
larisability by a factor (a/oo)^, where a is the typical size 
of the grain and aa = h? /{me^) is the Bohr radius. Later 
on, Rice 1 1 1 1 stressed the lack of screening in the calcu- 
lation of GE and showed that when screening is taken 
into account, one recovers the classical polarisability. We 
will see in the following that it is useful to consider the 
flux dependence of the polarisability within this crude 
approximation, since many results remain qualitatively 
true in the presence of screening. 



A. Position of the problem 

The polarisability of a system of non-interacting elec- 
trons can be understood as the sensitivity of the energy 
spectrum to an external electric field, just like persistent 
currents measure the sensitivity of the spectrum to an 
Aharonov-Bohm flux . As a matter of fact, the eigen- 
states |a > of the system and the eigenvalues Ca, are 
modified by an external electric field E. 



< a|eE.r|a > + X! 



< Q;|eE.r|/3 >| 
e,3 — Eq 



(14) 



E< (3\eE.r\a >,^ 
—-2 J 1/3 > 



/3#a 



e/3 - £0 



Consequently, it exists an induced charge density 
Sp{r) — e^^^^ IV'L('')I ^ lV'Q(r)| , associated with the 
asymmetry of each wave function generated by E. Using 
the definition ^ of a, the expression of the polarisability 
of a system of non-interacting electrons at zero temper- 
ature can be found: 



N 



2e' V V |<a|E.r|/?>| 
£;2 Z^ Z^ 



a=l,37^Q 



e/3 



(15) 



This expression depends on the eigen-energies of the un- 
perturbed system and its eigen-functions through the ma- 
trix elements of the position operator. Due to the energy 
denominator, the polarisability will in particular be very 
sensitive to the electric field induced coupling between 
the last occupied levels and the first non occupied ones. 
Note also that this expression of the polarisability is also 
from expression |l4| identical to the second derivative of 
the total energy of the system with respect to electric 
field. 



N 

E 



dE^ 



These results present strong similarities with the param- 
agnetic contribution of orbital susceptibility of a meso- 
scopic ring pierced by a flux line which is identical to ex- 
pression |l5| where the E.r is replaced by eA.p where A is 
the potential vector and P is the kinetic momentum op- 
erator. The total magnetic susceptibility contains also a 
constant diamagnetic term —Ne^/m ||l^, [Q which does 
not exist in the electric response. 



B. 1 dimensional Aharonov-Bohm ring 



Using expression |T^ it is possible to calculate exactly 
the induced dipolar moment in a Id non-disordered ring 
in an in plane electric field . 



Wavefunctions only depend on the angle , indicating 
the position in the ring, and satisfy the Id Shrodinger 
h^- a^Ma\ - T^Mo\ Furthermore, the 



equation: 



rV(6') = Eil){B). 



magnetic flux, associated with the periodicity in the ring, 
imposes the following boundary condition "0(0 + 27r) — 

e ^'^tp{9). This equation can be solved. Using the par- 
ity and periodicity of the spectrum, we then order the 
eigenvalues and the corresponding wavefunctions in as- 
cending order in the interval — ^ , ^ 



'=2p — 2mL2 
h 



'^2p+l 






^>'[~p+±^ 



(16) 



As a consequence of the particular geometry of the 
system, the electric field only couples certain states. In 
particular, the matrix element between adjacent states 
cancels unless the quantum numbers p et q are such that 
p — q = ±2, in this case : < p|X|(7 >= ^ The polaris- 
ability depends on the parity of the number of electrons 
N and reads : 




FIG. 2. Magnetopolarisability of a Id ring, of length 
L = 100, calculated with the Anderson model, for several 
values of the disorder: w = 10""; 0.1; 0.4; 0.7; 1.4. 



aw($) 



87r2 Vew+2(*) -ew($) 



eAr+l($) -eAr_i($) 



(17) 



At zero flux and for A^ ^ 1: aw — -^ — tt- It is 

o7r uq iv 

possible to evaluate S^a by computing (aw+i — aj^)/2. 
As a result : 
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(18) 



It is a very small effect, since it decreases rapidly with 
the number of electrons in the ring. This result remains 
true for a multichannel ring. We attribute this effect to 
selection rules inherent to the square lattice which are 
responsible for the cancellation of the matrix elements of 
operator X between eigenstates close to a level crossing. 
In particular it has been shown ||I3| that these selection 
rules do not exist in the hexagonal lattice where giant 
magnetopolarisability is expected for particular values of 
flux at the same level of approximation. 

In order to investigate the effect of disorder we have 
performed numerical simulations using the Anderson 
model. On figure g is plotted the magnetopolarisabil- 
ity for a ring of length L = 100 and several values of 
the disorder w. In absence of disorder, the results are in 
qualitative agreement with formula Eq. The effect then 
decreases when the disorder is increased, with a law close 
to :^. One is tempted to attribute this result to occur- 
rence of localization in the ID ring. However we will see 



in the following that this cancellation of the flux depen- 
dent polarisability is also observed in a multichannel ring 
or disk in the diffusive regime. 



C. Diffusive system 

In the following we discuss a diffusive system character- 
ized by a diffusion coefficient D — {l/d)vFle where le is 
the elastic mean free path which is assumed to be shorter 
than the system size a along the electric field. It has 
been shown 0, 111], [|l5| using semiclassical arguments 
that the average square matrix element, |< a|Ar|/3 >| 
depends mainly of the energy difference e = |e^ — eaj and 
the Thouless energy Ec = hD/L"^. For e < Ec it is of 
the order of a/g where g = Ec/A is the ratio between 
the Thouless energy and the mean level spacing A. At 
higher energy it decreases like 1/e^. This behavior is il- 
lustrated in fig.^ showing numerical results on a sample 
for different values of disorder where these two regimes 
can be clearly distinguished. 

It is then easy to deduce the order of magnitude of the 
polarisability noting that the summation in expression 
[l5| can be restricted to \a — P\ < g as a result: 



'a^A 



(19) 



The GE result is then recovered assuming that the 
Fermi wave length is of the order of the Bohr radius. We 
will see in the next section that this result is modified by 
screening. It is however worth continuing this analysis in 




FIG. 3. Energy dependence of the non diagonal matrix el- 
ements in a ring of length L = 60 with M = 4 channels, and 
several values of the disorderiu; = 0.9; 1.2; 1.4; 1.7. Note the 
additional step observed for e ~ 5A, which corresponds to 
the Thouless energy in the transverse direction of the ring: 
■^ where W is the width of the ring. 



this oversimplified picture of non interacting electrons in 
the case of an Aharonov-Bohm ring. From expression ll^ 
the polarisability is expected to exhibit flux dependence 
both from the matrix elements and energy denominators. 
We discuss in the following separately these two contri- 
butions which is justified in the context of random matrix 
theory where eigenfunctions and eigenenergies constitute 
two sets of independent random variables [|l6| . 
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1. Energy denominator 



We focus on the flux dependence of the quantity K — 



X]a</3 — - — )• This quantity depends only on the en- 
ergy differences (cq — e/3), and can be expressed as a func- 
tion of i?(e), the two levels correlation function : 



K 



de / du 



R{u) 



(20) 



In the diffusive regime, R{s) is well described by ran- 
dom matrix theory |]l6| and the average probability to 
find two degenerate adjacent levels is zero. This prop- 
erty is characteristic of level repulsion in the spectrum 
of a random matrix, which is stronger in a system where 
reversal symmetry is broken compared to a system where 
it is not. As a consequence, K decreases as a function of 
the flux at low magnetic flux and increases back in the 
vicinity of cf'o/'^ where time reversal symmetry is recov- 
ered. 

We have calculated numerically isT in a disordered ring 
of length L = 80 and width a = 8 for several value of the 
on-site disorder w corresponding to the diffusive regime. 
Figure |j illustrate the effect described above. In the ab- 
sence of other contributions, this term would give rise to 
a negative magnetopolarisability 



FIG. 5. Flux dependence of the diagonal and non-diagonal 
matrix elements in a ring of length L = 60 and M — 4 chan- 
nels, for w — 2. The factor 2 at zero flux between the diagonal 
and non-diagonal matrix elements is due to the fact that only 
the states such that ea < ep. 



2. Matrix element 

On figure ^ are plotted the fiux dependance of the av- 
erage square of the diagonal and non-diagonal matrix el- 
ements for a diffusive ring. It can be noted that they 
have opposite fiux dependences, resulting from the fact 
that Tr (X^) = Y.C. ^L + Ea#/3 Ki5 is ^ux indepen- 
dent. The diagonal elements decrease as a function of 
the flux, whereas the non-diagonal elements increases. 

The time reversal property of the operator X implies 
that its diagonal matrix elements are real even func- 
tion of flux and can be developed in successive pow- 
ers of cos(27r0/0o)- < ^aa > ^^ maximum for mul- 
tiple values of ^o/S. It is possible to evaluate analyt- 
ically this flux dependence of the diagonal matrix el- 
ements from random matrix theory, using the relation 
between Xa^a and the sensitivity of the energies to an 
electric field: Xa^a = {^e) f-o' S™"^^ electric field pre- 
serves time-reversal symmetry, the typical value of the 
derivative of the energy levels with respect to the elec- 




FIG. 6. Flux dependence of the matrix element and en- 
ergy denominators contributions to the average unscreened 
polarisability of a disordered square 30 x 30, calculated with 
the Anderson model for w — 3. These 2 contributions nearly 
compensate yielding to the absence of magnetopolarisability. 
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that this last result is qualitatively true 

1 



> is proportional to 4. Consequently: 

X,^ „ I (d) = 0) . Figure ^ shows 
Therefore, in 
absence of the energy denominator — ^ — , the flux depen- 
dence of the matrix elements would give rise to a positive 
magnetopolarisability. As a result this flux dependence 
is opposite to the contribution of the energy denomina- 
tors. Moreover we can see on fig. @ that these contri- 
butions almost exactly cancel and there is no magneto- 
polarisability in this model. 

We will see however in section V that a non zero mag- 
netopolarisability is found at finite frequency in the grand 
canonical Ensemble or at finite temperature in the canon- 
ical Ensemble where the energy denominators contribu- 
tion disappears. At this stage it is interesting to empha- 
size that the flux dependence of the average square of 
the X operator matrix element is opposite to the same 
quantity related to the current operator which changes 
sign by time reversal symmetry. So when computing the 
magnetic susceptibility the contribution of the matrix el- 
ements and energy denominators are of the same sign. 
This effect is related to the existence of a finite average 
current in the canonical ensemble for a diffusive ring. 

Note that both persistent current and zero frequency 
polarisability can be expressed as a function of the free 
energy of the system: / = —dF{E, 4>)/d4) and a ~ 
-d^F{E,(j))/dE'^. So: 



da d^L 



per 



dE^ 



(21) 



Let us emphasize that the absence of thermodynamic 
magnetopolarisability thus related to the insensitivity of 
persistent currents to an electrostatic field and can be 
qualitatively understood using a semi-classical argument. 
The effect of the electrostatic potential slowly varying at 
the scale of the Fermi wavelength can be included into 
the phase of the electronic wavef unctions. Persistent cur- 
rents can be computed from the integral of the classical 



action on a circular orbit. The contribution of the static 
electrostatic potential on such an orbit: {e/hvf) J V{r)ds 
where T^(r) = Ercos{9) and ds — rdd obviously cancels 
out. This argument is helpful to understand why the 
polarisability is independent of the flux in this very low 
frequency regime. 



IV. SCREENING 

A. The Thomas-Fermi approximation 

As mentioned above the mean field approximation re- 
duces the complicated many body problem of interacting 
electrons to a much simpler one, in which electrons move 
in an effective potential 0(r), resulting from the screening 
of the applied potential by other electrons. Therefore, it 
goes beyond a simple electrostatic calculation by taking 
into account the kinetic energy of electrons as well as 
their fermionic character [ [l7[ . 

In the linear regime, i.e for a small external potential, 
the induced charge density is simply proportional to the 
effective potential ||l|], fTgl: 

Spiv) = -e^n{EF)(Pir) = -eofc,''/'(r) 



at 3D with fc. 



sM^and 



Spir) = -e^n{EF)dp{v) = -e^k,cp{v)5{z) 

for a conducting plane at z=0 with ks = '^ ^ . 

The resolution of the self-consistent equation || within 
this Thomas Fermi approximation for a 3D sample in 
the presence of a uniform external applied field, gives 
rise to a screened potential whose value is of the order 
of E/ks which is confined to the border of the sample 
within Xg = 1/kg. This result is in principle not true 
any more in a 2D system where classical screening al- 
ready involves charge displacement in the whole system. 
However since the charge distribution is always singu- 
lar on the edge of the sample it is possible to approx- 
imate the screened potential for a disk of radius R by: 
<^{r,e) = {E/ks)F{{R - r))cos[e) where F is a peaked 
function centered on zero of width l/fc^. When screen- 
ing is taken into account the expression of the quantum 
polarisability at T=0 reads: 



2e^ 
£;2 



E 



< l3\(j>T.F{'r)\a X /3|E.r|Q; > 



(22) 



Within this approximation the GE result for the polar- 
isability is modified by a factor l/{aks)^ and becomes 
identical to the classical result. 



B. Beyond The Thomas-Fermi approximation 

As we have seen, the Thomas-Fermi approxim.ation 
does not take into account contribution of the screen- 
ing due to electronic interferences. They indeed give rise 
to a quantum correction for the response function x, and 
consequently to the effective potential and the polaris- 
ability a: 



X = XT.F + Sx 

(j) = <pT.F + 54> 
a = OLT.F + ^01 

where x is the response to the local field. Assuming 
that those corrections are small compared to the Thomas- 
Fermi value, it is possible to show, in agreement with 
Efetov [£2l (see appendix ^ for the detailed calculation) , 
that 5a can be simply expressed as a function of the 
mesoscopic correction to the one electron response func- 
tion 5x and the Thomas-Fermi potential: 



5a « -^Tr {(J)t.f5x4't.f) 



(23) 



In the following we discuss the response to a time depen- 
dent electric field. We will see that dynamical polaris- 
ability can be very different from the static one. 



V. EXPRESSION OF THE POLARISABILITY 



A. Response to time-dependent potential 



In addition, we assume that the effective potential is 
always in phase with the external one. This is justified 
if the dissipative part of the polarisability is very small 
compared to the non-dissipative part (a" <?C a'), this last 
quantity being discussed in |2l| . 

The expression of the density matrix is then obtained 
by solving the master equation. From it, we deduce the 
response function and later on the polarisability and its 
quantum correction 5a^, which can be expressed as a 
function of the eigenstates and eigenvalues of the hamil- 
tonian Hq. 




E 

Q#/3 



£« - e/3 £« - e/3 + ^ - «7 
7 v^ dfa 



— < P\4>T.F\a >< /3lE.r|a > 



"f + iuj 



y ^ < a\(t)T F\a >< a\E.r\a > ] (25) 

a / 



S^a 



2e^ 



S^Re 



E 



7 -|< l3\(f>T.F\a >\ 



7 -|- «w 



Y^^\<a\ct>T.FW>A (26) 

a " / 



The application of a time-dependent external poten- 
tial V{t) = eE.rexp(ia;t) raises the problem of the re- 
laxation of the system towards equilibrium. This process 
is made by inelastic processes characterized by a typical 
time scale, related to inelastic collisions. 

In the limit of a weak coupling, this process can be 
described by a master equation on the density matrix 
0: 



ih 



dp 

dt 



[Ho + V{t),p]-ij{p-p,g) 



(24) 



The parameter 7 represents the broadening of the en- 
ergy levels and characterizes dissipation processes, allow- 
ing the relaxation of the system towards equilibrium. The 
density matrix p^q satisfies the condition [H,peq] ~ 0, 
H = Hq + eE.r exp (iust) + (/)i„d(r, 7, t) being the hamil- 
tonian of the system in the mean field approximation, in 
an AC electric field at frequency uj. Two limits are to 
be considered. When the frequency of the time depen- 
dent potential is small compared to 7; the system follows 
the potential and stays at every moment at equilibrium. 
On the contrary, at high frequency compared to 7 , the 
system is always out of equilibrium. 



B. Quantum correction to the polarisabiHty in the 
diffusive regime for the different statistical ensemble 

In this section, we examine the differences between the 
canonical and grand canonical ensembles. The canonical 
ensemble corresponds to the situation for which the num- 
ber of electrons in the system is fixed. It is in particular 
the case for electrically isolated systems. The chemical 
potential is determined self-consistently by the condition 
N = J2ifo{^i — m)j and therefore depends on the flux 
through the energy levels. On the other hand, in the 
grand canonical ensemble, the system can exchange elec- 
trons with the thermodynamic reservoir, which imposes 
the value of the chemical potential. A physical realiza- 
tion consists in connecting the ring to a large metallic 
pad. At finite temperature, the occupation of the energy 
levels is spread over an energy interval of the order of 
the temperature. As a result, the sensitivity to the flux 
dependence of the chemical potential is suppressed when 
T ^ A, where A is the mean level spacing. Therefore 
the differences between CE and GCE are expected to dis- 
appear with increasing temperature and frequency. We 
first discuss the polarisability in the GCE. 




0.2 0.3 
*/*0 
FIG. 7. Dynamical {uj ^ 7)) and thermodynamical 
{u! <^ 7) polarisability in the grand canonical ensemble of 
a disordered ring 80 x 8, calculated with the Anderson model 
for w — 1 et w = 2. 



1. Grand canonical ensemble 

Due to the ergodicity property in the diffusive regime, 
the disorder average is equivalent to the average over the 
number of electrons. This property simplifies the calcu- 
lation. By averaging over the whole spectrum, one can 



show that < 



fc-ff) 



>= -i'^iicrc<>= ^/; 






Considering only regimes where w <C A and 7 <C A for- 
mula( Era) greatly simplifies: 



(5$ a 




7 



"f + iuj 



'T.Fia >\ 



]S{ea ~ y) |< a\(j)T.F\a >Y 



(27) 



of radius R and width W . The relaxation term cancels 
and only remains the term: 






GCE 




'T.F\a > 



(29) 



Therefore, one needs to evaluate the matrix element 
2 
. It is then possible to use the semiclassical cal- 



T.F 

culation of Mac Millan l2 



which yields to a general ex- 
pression of the average squares the matrix elements of 
the operator e"*'"" in a diffusive system which general- 
izes semi-classical expressions of the matrix elements of 
X relevant to the unscreened case discussed above. 



|<a|e''i-''|/3>f = 



Dq^ 



nh D'^q^ 



+ (eq - e/j) /h 



(30) 



In order to determine the relevant spacial frequencies 
of the screened potential, this latter is then decomposed 
in Fourier serie: 0t.f(i") = J2n=-Ai/2 ^ne"'"''' + e*'*"-'", 
where q^'" = ^Ur±^ug. At low energy, the main con- 
tribution to the matrix elements of the screened potential 
is dominated by the terms with the smallest wavevec- 
tor corresponding to n = 0, qmin 

A? ~ 1 _ A I 8R\s ) ^ 

""" hD 



~ ^ for which: 
;;^(l^)^> where A^ is the screen- 



ing length, Ec = 5;^ is the Thouless energy and A the 
mean level spacing between energy levels. 

We can then estimate the flux dependence of the ma- 
trix elements |0Q/3p as was done in section III for the 
unscreened potential: 

S^ < |0Q/3| >^= -Srj, < l^aal >^ 

= < |0qq| >GUE — < \4>aa\ >GOE 

Thus, the mesoscopic correction to the polarisability 
reads: 



Two frequency regimes can be distinguished, depending 
wether the frequency is smaller or larger than the relax- 
ation energy 7. In the limit lj ^ 7, which corresponds to 
a static electric field , the polarisability is simply equal to 
the flux dependence of the trace of the screened potential: 



Ep 



t>T.F) 



(28) 



Since the trace is independent of the basis which is con- 
sidered, this quantity is independent of the magnetic flux. 
Therefore, at low frequency, there is no mesoscopic cor- 
rection to the polarisability in the GCE as shown on flg- 
ure0. 

We now turn to the dynamical limit, ie for w 3> 7. As 
for the canonical ensemble, we consider a diffusive ring 



ao 



37r3 



Wg 



(32) 



where ao is the classical value of the polarisability Q g = 
-^ the dimensionless conductance. This result is in good 
agreement with the calculation of Blanter and Mirlin p^ ] 
using super-symmetry techiques. 



^We took the 2d limit corresponding to a disk for which 



oto 



cq-^R ■ For a quasi-ld ring, qq = 



[IB). Only the numerical factor in formula 
one consider this last value. 




(see section 
is modified if 



0.06 



0.04 - 
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isability reads in this ensemble: 



0.2 0.3 
*/*0 
FIG. 8. Dynamical polarisability (u) S> 7)) in the grand 
canonical ensemble in a disordered ring 80 x 8, calculated 
with the Anderson model for several value of the disorder: 
w = 1.; 1.5; 2; 3; 4. 



The magnetopolarisability increases when a magnetic 
flux is applied, corresponding to a positive magnetopo- 
larisability. A noticeable result is that this effect is 
inversely proportional to the conductance and thus in- 
creases with the disorder in the diffusive regime. It is 
expected however to decrease again with disorder in the 
localised regime when wave functions and eigen energies 
become insensitive to the Aharonov Bohm flux. 

This disorder dependence of — is illustrated on fig- 
ure g. On the other hand, the magnetopolarisability de- 
creases with the electron density rig, since the conduc- 
tance increases and the screening length decreases with 

The ring geometry is particularly favorable in order 
to observe this effect. The extrapolation of equation (|3^) 
leads to 5a/a ~ (A/2i?c)(As/a) for a two dimensional 
sample such as a disk or a square of typical size a and 
6a/a — {A/2Ec){Xs/R)'^ for a sphere. In a ring etched in 
a semiconductor heterojunction GaAs/GaAlAs with the 
following parameters: L = 8/xm, A^ = 400A, M — 10, 
Ec = 7A, we obtain 6a/a ~ 3 x 10"^. 



2. Canonical ensemble 

From relation ^, we deduce that at zero temperature 
and zero frequency the quantum correction to the polar- 



This last point is only true in 3d: The screening length 
decreases with Ue in 3d and is independent of n^ in 2d. 




< a\(t)T.Fir)\l3 >\ 

^13 — Cq 



(33) 



Where Nt is the total number of states. Note that this 
expression is very similar to the Gorkov Eliashberg one 
where matrix elements of the X operator are replaced 
by the matrix elements of the screened potential. Just 
as previously discussed in the absence of screening, the 
flux dependence of the energy denominators compensates 
exactly the one of the matrix element. However this com- 
pensation does not exist anymore, as it has been recently 
pointed out by Blanter et al. p4| when the frequency is 
of the order of the level spacing due to the disapearance 
of the contribution coming from the energy denomina- 
tors. A similar effect is obtained when increasing the 
temperature as shown precisely in R]. 



C. Effect of interactions in the canonical ensemble 

We have shown that in the canonical ensemble in the 
diffusive regime, there is no mesoscopic correction to the 
polarisability. However, we have neglected in this calcu- 
lation electron-electron interactions. 

In order to investigate a possible effect of electron- 
electron interaction a first and simple approach con- 
sists in taking the interaction potential as a perturba- 
tion. Rather than a Coulomb potential, we use an on- 
site interaction U5{y — r'). This approximation is justi- 
fied for high electronic density, since the interaction be- 
tween two electrons is then strongly screened by other 
electrons. In first order in U , the variation in total en- 
ergy can be expressed in terms of the local electronic 



density n(r) 



Z.^a=l 



< a r > 



SEtot = U / n{vYdr 



(34) 



The correction to the canonical polarisability is ob- 
tained then by the formula: 



Sur 



d\5Eto 
d^E 



(35) 



On figure ^ numerical calculations of this quantity 
Saint /oL is shown for a ring for various values of disor- 
der. 

Interactions give rise to a negative magnetopolarisabil- 
ity. The effect is more important in the ballistic regime 
than in the diffusive regime, in which it does not seem 
to depend on disorder. It gives rise to a magnetopolar- 
isability which is quite enhanced compared to its value 
in the absence of interactions. Nevertheless, the order of 
magnitude of this effect remains small compared to the 
dynamical polarisability. 
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A. Summary of the results 



s 




0,0 0,1 0,2 0,3 0,4 0,5 0,6 

FIG. 9. Correction to the canonical polarisability by inter- 
actions in a ring of length L = 60 and M — 4 channel, for 
several value of the disorder:^ = 0.9; 1.2; 1.4; 1.7. 



In the following table, are summarized the results ob- 
tained for — in the ballistic and diffusive regime, for the 
canonical and grand canonical ensemble. 



Ballistic regime 




t^ <7 


A >tj >7 


C 


-1 


-1 


GC 








Diffusive regime 




w ^ 7 


A >tj >7 


C 








GC 





C 8 ^ Af 1 

V Stt^ / W g 



VI. COMPARISON WITH THE EXPERIMENT 

The magnetopolarisability of two dimensional rings has 
been recently measured by Deblock et al |jl|. The sam- 
ple consists in an array of rings fabricated by electronic 
lithography in an heterojunction GaAs/GaAlAs. The au- 
thors use a resonant technique in which the rings are 
coupled to the capacitive part of a high frequency super- 
conducting resonator. 

The polarisability exhibits oscillations as a function 
of the flux, with a periodicity corresponding to (I)q/2 
through a ring. The order of magnitude as well as the 
sign of the effect is in good agreement with our grand 
canonical results. In addition, it decreases as a function 
of electronic density, according to the l/g dependence of 
formula ^ 

Taking into account the fact that the experiments 
where done at a frequency which is of the order of 1/3 
of the level spacing, where differences between canonical 
and grand canonical results are strongly reduced, this 
result is in good agreement with theoretical predictions 



VII. CONCLUSION 

We have calculated the average polarisability of meso- 
scopic rings and squares. In the grand canonical ensem- 
ble, we found that there exists a positive magnetopolaris- 
ability for frequencies larger than the typical broadening 
of the energy levels 7. The relative effect Sa/a scales like 
l/g, where g is the dimensionless conductance. In the 
canonical ensemble at zero frequency and zero temper- 
ature, the magnetopolarisability cancels in the diffusive, 
whereas in the ballistic regime a small negative effect is 
found. Differences between canonical and grand canon- 
ical ensemble disappear at frequencies or temperatures 
larger than the level spacing. Our results are in good 
agreement with recent experiments [fij. 



APPENDIX A: DETAILS ON THE 

CALCULATION OF THE MESOSCOPIC 

CORRECTION TO THE POLARISABILITY 

In this appendix, we derive in more details the ex- 
pression of the mesoscopic correction to the polarisabil- 
ity. In the following, we will note S = Stf + S'R and 
C = ^TF + S^, the response to the external and internal 
field respectively, defined in the matrix form as: 



Sp^X<t>^ 2(? 



(Al) 



The relation |l| which relate the induced charge density to 
the response function can be written in the matrix form: 



Sp = X(l - Ux) ^(t>ext 



(A2) 



The effective potential is also related to the response 
function by B[ The polarisability in the Thomas-Fermi 
approximation reads: 



Q-TF = ^Tr{xxTF<P) = ^Tr{xT.F{l - Uxt.fY 



E 



Tr{x:i, 



d<Pext) 



(A3) 



with: Ed — Xt.f{1 — Uxt.f)^^- The next step is to 
find the quantum correction ^S to the screened response 
function. We set by definition: 



E^x{l-Uxr'=^c 



-SE 



(A4) 



Assuming that quantum corrections are small compared 
to the Thomas-Fermi value, second order terms can be 
neglected: 



(1 - Uxtf)^ - USx^ci = xtf + Sx 



SE^Sxil^Ux)-\l + UEd) 
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Ep ^ 30K, our criteria allows to use linear response up 
to fields such that ER ~ 2mV. 



FIG. 10. Total energy of the electric field calculated in the 
Anderson model in a ring 40 x 5, and for a disorder w = 2. 

Using the fact that 1 + UEd ~ Xxp'^ch we show that 
the quantum correction to the response function can be 
expressed: 



6^ = ^ciXtfSxXtf^cI (A5) 

The quantum correction to the polarisability : 
6a ~ tr{x6'E,(j)ext) can then be written: 

Sa = tr{(j)TF5x4'TF) 

where we have used that XT_F'^ci0e2;t = 4>tf 



APPENDIX B: DISCUSSION OF THE LINEAR 
RESPONSE APPROXIMATION 

In order to estimate the validity of the linear response, 
we have calculated the total energy of a ring as a function 
of the electric field (screening is not here considered) . It 
can be expanded as: 



U{E) = Uo-doE-aE^ + ... 

The coefficient do correspond to the spontaneous dipo- 
lar moment of the ring, resulting from fluctuations of the 
charge in the presence of disorder. The ensemble average 
of do is zero. The energy average varies like E^ at small 
field. As it is shown on figure O , a deviation from the 
linear behaviour is hardly observed before a critical field 



eEjnaxR '^ Ef 



(Bl) 



where R is the radius of the ring. 

This critical value is less restrictive than the criteria B 
and justifies the use of the linear response up to very high 
fields. For instance, in a 2D electron gas, obtained in a 
semiconductor hetero junction GaAs/GaAlAs for which 
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